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When a particle moves around a Kerr black hole, it radiates gravitational waves. Some 
of these waves are absorbed by the black hole. We calculate such absorption of gravitational 
waves induced by a particle of mass in a circular orbit on an equatorial plane around a 
Kerr black hole of mass M. We assume that the velocity of the particle v is much smaller 
than the speed of light c and calculate the energy absorption rate analytically. We adopt an 
analytic technique for the Teukolsky equation developed by Mano, Suzuki and Takasugi. We 
obtain the energy absorption rate to 0{{v/cf') compared to the lowest order. We find that 
the black hole absorption occurs at 0{{v/c)^) beyond the Newtonian-quadrapole luminosity 
at infinity in the case when the black hole is rotating, which is 0({v/cf') lower than the 
non-rotating case. Using the energy absorption rate, we investigate its effects on the orbital 
evolution of coalescing compact binaries. 

§1. Introduction 

Among the possible sources of gravitational waves, coalescing compact binaries 
are the most promising candidates for detection by near-future, ground based laser 
interferometric detectors such as LIGO, VIRGO, GEO600 and TAMA. When a sig- 
nal of gravitational waves is detected, we will attempt to extract parameters for 
binaries, such as masses and spins etc., from inspiral wave forms using the matched 
filtering technique. In this method, parameters for binaries are determined by cross- 
correlating the noisy signal from detectors with theoretical templates. If the signal 
and the templates lose phase with each other by one cycle over ~ 10^ - 10^ cycles 
as the waves sweep through the interferometer band, their cross correlation will be 
significantly reduced. This means that, in order to extract information optimally, 
we need to make theoretical templates which are accurate to better than one cycle 
during an entire sweep through the interferometer's band. ^-^ 

The standard method to calculate inspiraling wave forms from coalescing bina- 
ries is the post-Newtonian expansion of the Einstein equations, in which the orbital 
velocity v of the binaries are assumed to be small compared to the speed of light. 
Since, for coalescing binaries, the orbital velocity is not so small when the frequency 
of gravitational waves is in LIGO/VIRGO band, it is necessary to carry the post- 
Newtonian expansion up to extremely high order va. v. A post-Newtonian wave 
generation formalism which can handle the high order calculation has been devel- 
oped by Blanchet, Damour and Iyer. Based on this formalism, calculations have 
been carried out up to post^/^-Newtonian order, or 0{v^) beyond the leading order 
quadrupole formula. ^-^ Another formalism has been developed up to 0{v^) by Will 
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and Wiseman. This formalism is based on that of Epstein and Wagoner. 

Another post-Newtonian expansion techniques based on black hole perturbation 
formalism has also been developed. In this analysis, one considers gravitational 
waves from a particle of mass fi orbiting a black hole of mass M when fi <^ M. 
Although this method is restricted to the case ji <C M, one can calculate very high 
order post-Newtonian corrections to gravitational waves using a relatively simple 
algorithm in contrast with the standard post-Newtonian analysis. The gravitational 
radiation at infinity was calculated using this method to 0{v^) by Tagoshi, Shibata, 
Tanaka and Sasaki in a Kerr black hole case, and to 0{v^^) by Tanaka, Tagoshi 
and Sasaki^) in a Schwarzschild black hole case. 

Despite these works, a study of gravitational waves which are absorbed by a 
black hole has been developed very little in the context of the post-Newtonian ap- 
proximation. The post-Newtonian approximation of the absorption of gravitational 
waves into the black hole horizon was first calculated by Poisson and Sasaki in the 
case when a test particle is in a circular orbit around a Schwarzschild black hole. In 
this case, since the effect of the black hole absorption appears at 0{v^), this effect is 
negligible for the orbital evolution of coalescing compact binaries in the above laser 
interferometer's band. However, as we shall see in this paper, black hole absorption 
appears at 0{v^) when a black hole is rotating. Thus it is important to investigate 
the effect of black hole absorption to the orbital evolution of coalescing compact 
binaries when at least one of the stars is a rotating black hole. 

In order to calculate the post-Newtonian expansion of ingoing gravitational 
waves into a Schwarzschild black hole, Poisson and Sasaki^) used two types of repre- 
sentations of a solution of the homogeneous Teukolsky equation. One is expressed in 
terms of spherical Bessel functions which can be used at large radius, and the other is 
expressed in terms of a hyper geometric function which can be used near the horizon. 
Then the two types of expressions are matched at some region where both formulas 
can be applied. They obtained formulas for a solution of the Teukolsky equation 
which can be used to calculate ingoing gravitational waves to 0{v^^), although they 
presented formulas for ingoing waves only to 0{v^). 

On the other hand, another analytic technique for the Teukolsky equation was 
found by Mano, Suzuki and Takasugi. ^'^^ Since this method is very powerful for the 
calculation of the post-Newtonian expansion of the Teukolsky equation, we adopt 
this method in this paper. 

This paper is organized as follows. In §2, we develop the Teukolsky equation. 
In §3, first we describe the analytic techniques for the Teukolsky equation and then 
we solve the Teukolsky equation analytically assuming the frequency of the wave is 
small. In §4, the energy absorption rate is calculated to 0{v^) compared to the lowest 
order. In §5, we consider the effect of black hole absorption on the orbital evolution 
of coalescing compact binaries. Section 6 is devoted to summary and discussion. 

Throughout this paper we use units in which c = G = 1. 
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§2. The Teukolsky formalism 

In the Teukolsky formalism, the gravitational perturbations of the Kerr black 
hole are described by a Newman-Penrose quantity ^4 = —Cais-ysn^fh^n^'ifi^. Here 
Cap-yS is the Weyl tensor, = ((r^+a^), -Z\, 0, a)/(2i:), = (ia sin 6*, 0, sin 6')/(V2(r+ 
ia cos 9)). We use the Boyer-Lindquist coordinates {t,r,6,(f)) and E = + cos^ 9, 
A = r^— 2Mr+a^, where a is the spin parameter of a Kerr black hole. We decompose 
^4 into Fourier-harmonic components according to 

(r - ia cos 0)V4 = E / dcoe-^^'+'^^'P -2Sem{0)RemUr). (2-1) 

The radial function i?£mu; and the angular function gSemi^) satisfy the Teukolsky 
equations with s = — 2 as 

^'i(Y-^)-V{r)Rem. = Te^., (2-2) 



dr \A dr 



1 d ( . ^ d ] 99.9^ (m — 2 cos 9y 
' sm 6'— > - a"^w"^ sm"^ ^ ^ ^ 



sin6'ci6'l c?(9 J sin^e* 

+ Aato cos ^ — 2 + 2mauj + A -2S£m = 0- (2'3) 

The potential is given by 

nr) = -^^±ifc^ + «.r + A, (2.4) 

where K = (r^ + a^)a; — ma and A is the eigenvalue of -28^^. The angular function 
sSim{d) is the spin- weighted spheroidal harmonic which can be normalized as 



/ I^S^^p sin 6'd6' = 1. (2-5) 
Jo 



The source term T^^^^ is given in a paper TSTS. 

We define two kinds of homogeneous solutions of the radial Teukolsky equation: 

nin _^ f Cimuj^'^e-'''''* for r ^ r+ „x 

f R°"t pifcr* I /\2 pin -ikr* f 
^""^ iL'wrV'^^* for +00, ^ ^ 

where k = lu — ma/2Mr_|_, and r* is the tortoise coordinate defined by 

dr* , 
— — ar 
dr 

2Mr+ , r - r+ 2Mr_ ^ r - r- 

= '^+ —In In— r— , (2-8) 

r+ - r_ 2M r+ - r_ 2M ' ^ ^ 
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where r± = M ± \/ M'^ — a^. 

We solve the Teukolsky equation by using the Green function method. A solution 
of the Teukolsky equation which has a purely outgoing property at infinity and a 
purely ingoing property at the horizon is given by 

Rim^ = ^ |i?Zu, £ dT'Ftl^T,^^A-^ + /~ dr'RZjr,^^A-^^ , (2-9) 

where the Wronskian W is given by 

W = 2iiuDe„^coA'l^nu.. (2-10) 
Then, ip^ has an asymptotic property at the horizon such that 

RimUr ^ r+) ^ ZZ^A^e-"''*, (2-11) 

where 

= 2,J'^^^A^'^ r dr'RZ.T,m.^-'- (2-12) 
The energy flux formula is given by Teukolsky and Press as 
dE,,,,,\ _^ f,^ 2SL 128uk{P + 4€^)(fc^ + 16e^)(2Mr+)5 , , 

1^ \ \^lm^\ > l2-icij 



dtdn ) 27r |C| 

where e = /«/(4r+) and 

|Cp = (Q2 + 4aa;m _ 4a2a;2)[(Q _ 2)^ + 36aa;m - Sea^a;^] 



+(2Q - l)(96a^a;^ - 48aa;m) + 144a;^(M^ - a^), (2-14) 

and Q = A + 2. 

The post-Newtonian expansion of these formulas can be obtained by solving the 
Teukolsky equation by setting e = 2Muj ^ 0{v'^), z = tor ^ 0{v) and by assuming 
e ^ z ^ 1. Here, we define v = (M/ro)^/^, where ro is the orbital radius of the test 
particle. 

The solution of the spin weighted spheroidal harmonics was given explicitly in a 
previous paper up to O(e^) which is accurate enough for the present purpose. In 
the next section, we explain the method to obtain the homogeneous solutions of the 
radial equation. 



§3. Analytic solutions of the homogeneous Teukolsky equation and the 

post-Newtonian expansion 

New analytic representations of the homogeneous Teukolsky equation were found 
by Mano, Suzuki and Takasugi. ^'^^ ' In this paper, we adopt variables used by Mano 
and Takasugi. Since all the details are given in that paper, we only describe the 
outline of the method here. In this section, the value of the spin s is assumed to be 
s = — 2. In this method, the solutions of the Teukolsky equation are represented by 
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two kinds of expansions. One is a series of hypergeometric functions which is defined 
as 

^ r{l-s- 2ie+)r{2n + 2v + l) ^ 
^ r{n + v + l-iT)r{n + v + 1- s-ie)"''^ 

n=—oo \ ' • / \ • • J 

xx"'F ^— n — u — IT, —n — u — s — ie, —2n — 2v; 4^ , (3-1) 
where F{a, b, c, x) is the hypergeometric function, and 



The coefficients obey the three terms recurrence relation 

«+l + /?>n + Tn^n-l = 0> (3-3) 



where 



ieK{n + v + 1 + s + ie){n + v + 1 + s — ie){n + v + 1 + ir) 



(n + i^ + l)(2n + 2zy + 3) 
= -A - s{s + 1) + (n + + + 1) + + e(e - mq) 
e(e — mq){s'^ + e^) 



ien^n + v — s + ie){n + u — s — ie){n + v — ir) fo a\ 

" (n + j/)(2n + 2:/- 1) ' ^ ' 

The series converges if v satisfies the equation 

i?„(i/)L„_i(i/) = 1, (3-5) 
where Rn{i^) and Ln{i') are continued fractions defined by 

= = "n-l7n an-27n+l __ _ (g.g^ 

Then the expansion equation (3-1) converges if r < oo. 
The ingoing solution R^^^^ is given by 

-^w = e*'"'(i?o + ^o"'^)- (3-7) 
The other solution is a series of Coulomb wave functions, defined as 
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where (a)„ = r{n + a)/r{a), and Fn+us{z) is the Coulomb wave function defined 

by 



^ ^ ' r(2n + 2z/ + 2) 

x$(n + z/ + 1 - s + ze, 2n + 2z/ + 2; 2?2;), 

where $(a, 6, 2) is a confluent hypergeometric function, z = u{r — r_) and e± = 
(e lb t)/2. The expansion coefficients a':^ obey the same recurrence relation (3-3). 
The expansion equation (3-8) converges if is a solution of Eq. (3-5). In that case, 
the series converge for z > en. 

This solution is decomposed into a pair of solutions, the incoming solution at 
infinity, -R+;s, and the outgoing solution at infinity, RL-^g- Explicitly, we have 



R'c;s=K;s + R%s. (3-9) 



where 



r{y + 1 + s-ie) ^ ' 

00 

X ^ i"<(22)"^'(n + z^ + l-s + ze,2n + 2i/ + 2;2i5), (3-10) 

n=— 00 

00 
n=— 00 

X [''^|^^~'^|" a;;(2gf^(n + z^ + 1 + s - 2n + 2z/ + 2; -2?5),(3-ll) 

where ^{a,b,z) is a confluent hypergeometric function. The upgoing solution 
is then given by 

I, 

Hmu) 

We see that the above two kinds of solutions have a very wide convergence region. 
We also see that the two solutions Rq.^ and Rq.^ are not linearly independent of each 
other. Then Rq.^ must be proportional to R(j.g- The relation between them is given 

by ' 

where 



{2eKy+'-^2-H^r{l -s- 2i€+)r{f + 2z^ + l)r(f + 2z/ + 2) 

r{iy + 1 - iT)r{iy + 1 - s - ie)r{f + u + 1 - s + ie) 
r{p + 1 + iT)r{u +1 + s + ie) 



X 



r(f + U + 1 + iT)r{r + 1^ + 1 + s + ie) 
f + 2z/ + l)„(i/ + l + s-k)„ ^ 



Coo 
(n-r)! (z/ + l-s + ze)„ 

xf E M^.^'^o o^ ^''^]^'"'\^<is)y\ (3-13) 

V^fr'oo - + 2z. + 2)„ (z. + 1 - s + ze)„ ' ^ ^ 
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where f can be any integer and K^i^s) is independent of the choice of f. This relation 
holds in a region where both expansion equations (3-1) and (3-8) converge. 

Using this relation, it is possible to obtain the asymptotic amplitudes of 
at infinity. They are given by 



A out 



LU 



K,{s) 



le 



, sin7r(i^ — s + ie) 
simriv + s — ie) 



v{s) 



l+2s 



K,{s) + ie'''''K^,-i{s) 



ie In e 



where 



2-l+s-iegi(7r/2)(i^+l-s)g-7re/2 



-1 — s+ie 



-j(7r/2)(i/+l+s) -7re/2 



r{i^ + l- s + ie) 
r{i/ + l + s-ie) 

{u + l + s 



oo 

E ' 

n=—oo 



E (-1) 



_ {i^ + 1 - s + ie)n 
The asymptotic amplitude at the outer horizon C^.^ is given by 



The asymptotic amplitude 



2s 



gie+ In K 



E 



of at infinity is also given by 



oo 

E 



(u + l + s 



(u + 1 — s + ie 

— no ^ 



,(3-14) 
(3-15) 

(3-16) 
(3-17) 

(3-18) 

In e 



UJ 



l+2s ■ 

(3-19) 



Next, we explain the relation of the method above to the post-Newtonian ex- 
pansion. It is easy to see that if we set Oq = 1 at n = 0, the order of in e increases 
with \n\. Basically, this is because a'^^ and 7^ contain an overall factor e. Further, 
from the behavior of each term of the expansion equation (3-8) as a function of z, we 
see that the post-Newtonian order of the expansion equation (3-8) increases. This 
can be checked by setting e ~ f and z ~ f , where v = (M/ro)^/^, and ro is the 
orbital radius of the test particle. Then, the above expansion is very useful for the 
post-Newtonian expansion. Note that these increases of the post-Newtonian order 
are not monotonic when n < 0. Therefore we must be careful to treat the expan- 
sion for n < 0. This is not, however, a serious problem because we can obtain the 
desired post-Newtonian accuracy by summing up the terms up to the appropriately 
large \n\. 

Finally, we briefly explain the procedure of the analytic calculation. We first 
calculate the eigenvalue A in a power series of e using a method due to Fackerell 
and Grossman. Next we determine u by solving Eq. (3-5). This equation can 
also be solved in a power series of e. Then we evaluate the expansion coefficients 
a'^. This can be done conveniently by using continued fractions Eqs. (3-6) and 
setting = 1. Once we obtain the a'^, it is straightforward to calculate asymptotic 
amplitudes A'^^, Cimui and D^,^. It is also straightforward but tedious to calculate 
from Eq. (3-11) by assuming e ^ 2; ^ 1. These results are given in Appendix 

B. 
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§4. The black hole absorption to 0{v^) 

In this section, wc evaluate Eq. (2-13). The calculation is almost the same as 
those given in section III of a paper TSTS. First, we solve the geodesic equation for 
circular motion in the equatorial plane. Next, we calculate the source term Timuj. 
The result is given by 



'-imuj 



f 



(^n n + n 



+ {{Afnnl + Afnml)Sir - To)} ^ + {Afnfh2S{r - Vq)} . 



e=ir/2 



, (4-1) 



where AnnO, etc., are given in Appendix A. Inserting Eq. (4-1) into Eq. (2-12), we 



obtain Zf^ as 



tniio f /I \ A 



dr 



'-{A 



mnl 



mml 



j2 pup 



dr"^ 



mm 2 



r=ro,0=7r/2 



= 5{u - mQ)Zf^ , 
where i? is the orbital angular frequency given by 

r 



Q = 



3/2 



(4-2) 



(4-3) 



Prom Eqs. (2-13) and (4-2), the time averaged energy absorption rate becomes 



'dtJii 



E 

em 

E 



(-) 

\dt L 



(4-4) 



Using a property of -25""^ (0) at ^ = 7r/2, it is straightforward to show that Te^^m-ui 
= {—lYT£^m,Lo where Tg^m,Lo is the complex conjugate of T^^m.^;- Since the homo- 
geneous Teukolsky equation is invariant under complex conjugation followed by 
m — > —m and u> — ^ —to, we have _^ _^ = (— Then, from Eq. (4-4), we 
have {dE/dt)i^^rn = {dE/dt)e,m- 

In order to express the post-Newtonian corrections to the black hole absorption, 
we define 77^^ as 



fdE\ _ 1 (dE\ 5 H 



\dt) 



dt J 



N 



e,m 2 

where {dE/dt)^ is the Newtonian quadrupole luminosity at infinity, 

_ OZ,/LX 1V± _ / ^ 1 ^ 

\dtj, ~ ' ~~ ~ 



(4-5) 



N 



M 



,10 



Post-Newtonian Expansion of Gravitational Waves 
In Appendix C, we give rjim for m > 0. The results for m < are given by rji^-, 

m-m- 

The total absorption rate is given by 



fdE\ JdE\ 
Kdt Ju \dt J N 



3 Q 1 / 33 o \ n 



(7 1 85 1 13 \ 

^q'' + 2qB2 + -+6q^B2 + —q^ + Sq\+-K+ — Kq^j 

/ 4651 o 43 17 A 4 

/569 2 371 4 o 3 o „ „ 33 4 

+ q^ + —q^ + 18 q^B2 - - q^Bi +2K + 2 + — q^K 

+6qB2 + ^Kq'^ + qBi^v^ 

f 2718629 , „ 428 2 ^ 428 , „ , ^ 
+ q-4B2 H 70 + -TT^'q-l (7ln2-4g(7s 

V 44100 ^ ^ 105 3 ^105^ ^ 

428 428 

-12 g3 C2 - 36 g^Ss - 56 g^Sa + — + — g3 In 2 + 2 q^w'^ 

35 35 

428 , 428 , 428 o , ^ 428 , , „ „ o 
+T05^^^'^ + T05^^^ + ^^ ^^^^ + ^7 + ^^ ^^-'^^^ 

856 , 856 o , B2 „^ g3 „ g 

-24 g^^a^ + g In + q^\iiv-4—-32—-31- 

105 35 K K K 

+57 — + g^/c + 5^5^ - ^Kq-^Kq^ -^q'^Bi 

0^52 g'^^g q^C2 qC2 q^B2 
-48 ^ — - + 28 ^ — - - 24 ^ — - - 8 — + 24 - — - 

Kj Ki tv tv tv 

2400247 o 299 A . 

q-^ H q' 

19600 ^ 16 V 

/225 . 41 R 86 8741 . 3485 4 167 « 
+ g 5? Kq*" -\ \ q^ H g^ H g^ 

V 28 ^ ^ 28 ^ 7 56 ^ 42 ^ 112^ 



-1 
+ (- 



86 45 9 899 803 o 1665 , 
— K-\ qBs q^Bi H qBi q'^Bi H q'^Bs 

7 56 ^ ^ 28 ^ ^ 168 ^ 224 ^ ^ 224 ^ 

2372 3 „ 16 5 „ 719 4 22201 , 796 „ \ 7 

q'^B2 q B2 H q k-{ Kq^ H gSa ]v' 

2j y ^ J y ^ ^ 12 ^ 168 ^ 21 ^ 

/ 20542807 „ 1061 13 , 995 , „ 

q - 12 B2 -2 Bi -\ 7 g H Tr^g H g In 2 

V 88200 ^ 35 '^ 6 ^ 21^ 

,^ ^ 308 4 „ 1496 2„ 12197 , 3873 „ 

-12 gC^ - 36 g^ - q^B2 - q^Bs + -j^ q^^ + ^ g^ In 2 

47 , 2 1391 , 428 , 5029 o , 37 g^ 1284 . ^ 
+ ^ rvr^ + --^qlnK + — qAz + — - g^lnK + — ^ + — — gM^ 

8 105 35 140 6 k 35 

9 00 4574 8613 o r„ q^i ^3 

-24g5,2 - 72g35,2 + __ ^In^; + ^ g^ In^; - 12 ■ 



341 g^ 

K 4 K 
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637 741 g5 73 283 . 

7- - + -j \- -^q Bi - qCi - —q Bi 

6 K 4 K 6 18 



-3^2 107 



-2 



+84 



3 

2" ' 105 
, 3 q^Ci 



qAj 



107 



140 



^q^Bi ^ 13 q^Bi 3 q^Bi 



+ 

K 

q*B2 



Bi 



2qBi^ + -q^Ci-2 
2 k 4 K 

q^C2 qC2 q^B2 



144 



q'B, 



2 



K 



K 



K 



2945984497 
6350400 



+- 



1385 



25 



(4-6) 



where 



B„ = 



1 
2 
1 

2i 



(0) 



mq 



Cf) — — 
2 



V^(i) (^3 + 



mq 



vT 



(^3 

+ f3- 



mq 



niq 



niq 



(4-7) 



and tp^'"\z) is the polygamma function. 

We see that the absorption effect starts at 0{v^) beyond the quadrapole formula 
in the case q 0. If we set q = 0, the above formula reduces to 



dE\ fdE\ f o 



10 



(4-, 



which was derived by Poisson and Sasaki. Then, the black hole absorption is more 
important in the case g 7^ 0. Wc also notice that Eq. (4-6) can be negative when 
q> Q. This is an example of super radiance phenomenon. 

It is not clear that Eq. (4-6) has a finite limit as \q\ — 1. However by using the 
formulas 



hm 3 + 



mc 



V^-q 



= ] =lnn-lnAv + iy^|, 



we can obtain the limit of {dE/dt)-^ as 



(4-9) 

(4-10) 



lim 



dE\ 



dE 



9-*±i \ dt Jji \ dt J ]\[ 



q 49 g 2 133\ 3 



6817 q 4 /535 97 ^ 5 

-r-rv + { \ Tr]v^ 

336 \q\ V 16 8 ' 
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/3424 q , , , 1712 q 3424 q , , , 

_ 3647533 289 , 

22050 a fi « ' / 



22050 \q\ 6 |g| 

/ 84955 55873 \ 7 

+ \ vr f ^ 

V 336 672 J 



/14077 q 16441 q 34987 g 193 2 9 

_4057965601 ^ _ 1289 ^ \ J ^ 
6350400 \q\ 9 \q\ J ' ^ ^ 



In Appendix D, we give {dE/dt)t{ written in terms of x = (Mil)^/^. 

§5. The orbital evolution of coalescing compact binaries 

Using the above results, we estimate the effect of black hole absorption on the 
orbital evolution of the coalescing compact binaries. We ignore the finite mass effect 
in the post-Newtonian formulas and interpret M as the total mass and /x as the 
reduced mass of the system. 

The total cycle N of gravitational waves from an inspiraling binary is given by 

^ QdE/dv 

N= I dv T7TT-r: (5-1) 

Jvs TT dE dt ' ^ ' 



where vi = (Af/r.;)-'^/^, Vf = {M /r fYl'^ . and r,j and Vf arc the initial and final orbital 
separation of the binary. The energy loss rate is given by dE/dt = {dE/dt)\i+ 
{dE/dt)oo, where {dE/dt)oo is the luminosity at infinity. 

We set Tf = 6M and define to be the radius at which the frequency of the 
wave is lOHz. We evaluate AN which is the difference of the total cycle N caused 
by the formulas for dE/dt with and without the black hole absorption effect. The 
value of AN is calculated explicitly by 



AN 



n dE/dv _ n dE/dv 

Jvf ^ TT (dE/dt)n + idE/dt)^ Jy, ^' 



vf TT {dE/dt)B. + {dE/dt)oo Jv; tt {dE/dt)^ 



(5-2) 



The results for the typical NS-BH systems are given in Table I. We find that 
the effect of black hole absorption is very small when g = 0. We also find that the 
black hole absorption is more important in cases for which q > than in cases for 
which q < Q. This is because, when g < the total cycle is much smaller than the 
case 5 > 0. We also notice that the black hole absorption is more important when 
the mass of the black hole is large. The reason is as follows. ^) 

When the black hole mass becomes larger, the initial non-dimensional orbital 
radius ri/M becomes smaller and thus Vi becomes larger. When g > 0, this effect 
is dominant. Hence, although the total cycle N becomes small when the black hole 
mass is large, the effect of the black hole absorption becomes larger. 
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Table I. The difference of the total cycle N caused by the formulas for dE/dt with and without 
the black hole absorption effect for typical neutron star(NS)-black hole(BH) binaries with mass 
(Mns,Mbh). 



q 


(1.4M0,1OM0) 


(1.4Mq,40Mq) 


(1.4Mq,70Mq) 





< 10-^ 


0.02 


0.1 


0.9 


0.5 


1.1 


1.5 


-0.9 


0.3 


0.5 


0.4 



§6. Summary 

We have calculated expressions for the gravitational waves which are absorbed 
by a rotating black hole and obtained the energy absorption rate by the black hole 
to 0{v^^) beyond the Newtonian quadrapole luminosity at infinity. We adopted a 
method developed by Mano et al. to calculate the post-Newtonian expansion of 
the Teukolsky equation analytically. We find that black hole absorption occurs at 
0{v^) beyond the quadrapole formula. This order is 0{v^) lower than in the non- 
rotating case. The energy absorption rate obtained here displays a property of super 
radiance. 

Using the formulas obtained above, we have estimated the effect of the black 
hole absorption on the orbital evolution of coalescing neutron star-black hole bina- 
ries. We calculated its effect on the cycle of waves from the binaries in the laser 
interferometer's band. We found that the effect of the black hole absorption is more 
important in the case when g > 0. (This is the case when a particle moves in the 
direction of the black hole rotation.) We also found that black hole absorption is 
important in the case when the mass of black hole is large. 

Gal'tsov ^^-^ gave formulas for black hole absorption from the radiation reaction 
force on the particle. Gal'tsov's formula is expressed as 




This formula agrees with our formula Eq. (4-6) to lowest order. 

In this paper, we have calculated the post-Newtonian formulas only to 0{v^) 
beyond the lowest order. However, it is straightforward to calculate higher orders 
with our method. 

The analysis in this paper has been restricted to the case in which a test particle 
moves in a circular orbit on the equatorial plane. However, it has been suggested 
that the inclination of the orbital plane from the equatorial plane will significantly 
affect the orbital phase evolution. Hence, the effect of the orbital inclination should 
be taken into account in future works. 
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Appendix A 

The Functions in 



In this appendix, we give functions which appear in Z^^. 
The functions A in Eq. (4-1) are given by 



^nnO 



Cnnp-^p-^Lt{p-^L^{p'S)}, 



27rZ\2 
2 



^/^^A 



CmnP 



(Lts) (^ + P + P 



—a sin 9S—{p — p) 



A 



mmO 



P ^pCfhmS 



-i 



A, 



mnl 



■^mm 1 — 



27r 
1 



p ^CfnnlLtS -\-iasme{p- p)S], 



P '"pCrhmSli^ + p 



A 



P pCmmS, 



where S denotes -28^, 

The functions C in the above formulas are given by 



P 



r - 



l2 



al 



Cn 



PP 

2Ei 



i sin 9 aE 



I. 



sm 



i sin 9 1 aE — 



sm^9 



(A-l) 



and i is given by the geodesic equations of the particle as 

r^i = -(aE- a + (^E{r^ + a^) - aZ,) . (A-2) 

The specific energy E and the angular momentum are given by 

1 — 2^2 + qv^ 
" (l-3^;2 + 2gt>3)i/2 ' 
_ rpyjl - 2qv^ + g^z;^) 

(1 - 3^;2 + 2gt,3)i/2 • ^"^'^^ 
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Appendix B 

The Functions m S 5 ■ 
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960 8 ^ ^ 112 ^ V ^ K ; 8 
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'4 189 ^ 49 3 3969 



1 1 67 , 
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Appendix C 

The Energy Absorption Rate for Each Mode 

The energy absorption rate for each mode to 0{v^) are given by 
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Appendix D 

The Alternative Form of {dE/dt)}i 



In this section, we describe the absorption rate {dE/dt)i{ by means oi x = 
(M!7)V3. Here we define 
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